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Abstract 

A class of exact general relativistic rotating thick disks surrounded 
by rotating dust is constructed by applying the "displace, cut, fill and 
reflect" method on the van Stockum class of metrics. Two particular 
solutions are considered in detail. We find that the disks have annular 
regions with negative energy density and heat flow is present in most 
of the radial extension of the disk. It is shown that thin disks are 
composed of exotic matter with an equation of state of cosmic strings 
or struts. We also comment the recent relativistic galactic model pro- 
posed by Cooperstock and Tieu and the presence of an additional disk 
of exotic matter in their model. 
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1 Introduction 

Exact solutions of Einstein's field equations describing stationary or static 
axially symmetric configurations of matter are of great astrophysical inter- 
est, since the natural shape of an isolated self- gravitating fluid is axially 
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symmetric. In particular, disk-like solutions can be used as models of galax- 
ies and accretion disks. Solutions for static thin disks without radial pressure 
were first studied by Bonnor and Sackfield [T] and Morgan and Morgan [2], 
and with radial pressure by Morgan and Morgan [3j . Several classes of exact 
solutions of the Einstein field equations corresponding to static thin disks 
with or without radial pressure have been obtained by different authors [Ij— 
[13) . Thin rotating disks that can be considered as a source of the Kerr 
metric were presented by Bicak and Ledvinka [T3] , while rotating disks with 
heat flow were studied by Gonzalez and Letelier [15]. Also thin disks with 
radial tension [16], magnetic fields [T7] and magnetic and electric fields [18j 
were considered. The nonlinear superposition of a disk and a black hole was 
first obtained by Lemos and Letelier |7j. Perfect fiuid disks with halos [19j 
and charged perfect fluid disks [20] were studied by Vogt and Letelier. The 
stability of some General Relativistic thin disk models using a first order 
perturbation of the energy- momentum tensor was investigated by Ujevic 
and Letelier [21]. For a survey on self gravitating relativistic disks, see for 
instance [22j . 

Even though in a first approximation thin disks can be used as useful 
models, in more realistic models the thickness of the disk should be consid- 
ered. The addition of a new dimension may change the dynamical properties 
of the disk source, e. g., its stability. Thick static relativistic disks in various 
coordinate systems were presented by Gonzalez and Letelier [23] and Vogt 
and Letelier [2l]. General Relativistic versions of some well known Newto- 
nian three-dimensional galactic potential-density pairs were studied by Vogt 
and Letelier [25j . 

In the works cited above an inverse style method was used to solve the 
Einstein equations, i. e., the energy-momentum tensor is computed from 
the metric representing the disk. Another approach to generate disks is by 
solving the Einstein equations given a source (energy- momentum tensor). 
Essentialy, they are obtained by solving a Riemann-Hilbert problem and are 
highly nontrivial. This has been used by a German group to generate several 
exact solutions of disks [26]-[33]. 

Recently Cooperstock and Tieu [34] devised a new approach to modelling 
galactic dynamics via General Relativity. A galaxy is idealised as an axi- 
ally symmetric uniform rotating fluid without pressure (dust). The authors 
show that the density distribution of the galaxy is related to the rotation 
parameter in a nonlinear way, and then show that several measured rotation 
curves for galaxies are consistent with the mass distribution of matter in flat- 
tened disks, therefore excluding the need for a massive halo of dark matter. 
However, their model has been criticized by some authors [35^ [36| I37j . In 



2 



particular, Korzyiiski [35j has argued that the model possesses an additional 
singular rotating thin disk at z = and Vogt and Letelier [36] have shown 
that this singular disk is made of exotic matter either like cosmic strings 
or struts with negative energy density. Balasin and Grumiller [38], using a 
somewhat different mathematical formulation, have shown that is possible 
to obtain flat rotation curves without unphysical sources at z = 0, but in 
turn distributional sources may arise on the z-axis. Motivated by Cooper- 
stock and Tieu's work, we investigate the properties of thick rotating disks 
immersed in rotating dust obtained by using the inverse style method (image 
method). 

The work is divided as follows. In Sec.[2]we present the van Stockum [39] 
class of metrics that will be used to construct the disks and the formalism 
to calculate their physical variables. A class of rotating disks based on the 
Bonnor [iQ] solution is presented in Sec.[3l In Sec. U] another class of rotating 
disks based on an extension of Bonnor's solution is discussed. The thin disk 
limit and some comments about the Cooperstock and Tieu galactic model 
are presented in Sec. [H Sec. [6] is devoted to the discussion of results. 



2 The van Stockum class of metrics and construc- 
tion of disks 

The class of solutions found by van Stockum ^39j describing stationary axi- 
ally symmetric distributions of dust matter has a line element given by 

ds2 = {dt - Ndiff - r^d^^ - e-'idr^ + dz^), (1) 

where N and v are functions of the quasi-cylindrical coordinates r and z. 
For metric Eq. ([T|), Einstein's field equations Gf^^ = {S-kG / c^)T^y can be 
cast as 

N,rr + Nzz - ^ = 0, (2a) 
r 



(2b) 



where p denotes the energy density of dust. Eq. (i2a|) can be expressed as 

<^,rr + + ^ = 0, (3) 

r 
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by the transformation N = r^^r- 

In order to obtain a thick disk surrounded by dust, given a solution of 
Eqs. ()2ap - (pc]) . we apply the transformation z h{z) + b, where 6 is a 
positive constant and h{z) an even function of z such that h{0) = 0. The 
parameter b controls how much matter is concentrated near the symmetry 
axis (larger values of b generate smoother distributions of matter); whereas 
the function h{z) dictates the matter distribution along the z direction. 
This is equivalent to construct a disk by the "displace, cut, fill and reflect" 



method that is described at length in 



The nonzero components of the 



energy-momentum tensor for metric equation ([T|) are 



levrGe'' 
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N.rr + N,, 
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(4b) 
(4c) 

ZZ ) 

(4d) 
(4e) 

(4f) 



Applying the transformation z h{z) + b and using Eqs. (j2ap - ([2cj) we ob- 
tain the following expressions for the components of the energy-momentum 
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tensor of the disk 



= poc^ + 



h I — ^ - 1 + (1 - /i ) I ly^hh - 
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(5b) 
(5c) 



+ Z^.h - (1 - /l ) o + 



2r2 



2 /L'S 



327rGr2e 



1) 



0, 



(5d) 
(5e) 



where primes indicate differentiation with respect to z and po is given by 
Eq. ^ with replaced by N^h- 

The physical variables of the disk are obtained by solving the eigenvalue 
problem T^^^ = A^", and has the solutions 



x± = -{s±Vd), 

S = Tl + = poc' + 



-h"v^h + {l-h 



12 \ 



2r2 



Poc 



Xr = , Xz = T. 



(6) 
(7) 



(8) 
(9) 



Defining an orthonormal tetrad (!/", VF", X", F") with V = 7Vi(l, T, 0, 0), 
W = A^2(A, 1,0,0), = 6-^^/2^0,0,1,0) and = e^^/^^o, o, 0, 1), where 
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A^i , are normalization factors and 



^ l(A,-r/)/ri D>o 

\{Tl-T^)/{2T^). D<0 

[0, D<0 
the energy-momentum tensor can be written in the canonical form 

Tab = pVaVb + P^WaWb + K{VaWb + WaVb) + Pr-XaXb + P.YaYb. (12) 

The energy density p, the azimuthal stress P^, the heat flow function k, the 
radial pressure Pr and the tension Pz along z are, respectively, 

(13) 




(15) 

Pr = -r;, p, = t;. (le) 

For the functions h{z) we use a class of even polynomials derived in [2 

h{z): 



—z + C, z < —a, 

Az^ + 5z2«+2^ -a<z<a, (17) 
z + C, z > a, 



with 



^ 2n + 1 — ad ^ ad — 1 ^ a{2n + 1 + ad) 

~ Ana ' ~ 4n(n + l)a2"+i ' ~ 4(n + 1) 

Here n = 1,2,...; a is the half-thickness of the disk and d is the jump of the 
second derivative on z = ±a. We have |/i'(2)| = 1 and h" = for \z\ > a, 
so Eqs. ([5a|) - p^ are reduced to van Stockum's solution outside the disk. 
For l^l < a we always have |/i'(-z)| < 1 and Eqs. (j5e|) and ()16p show that we 
have radial pressures and tensions along z. Also when d = h" {\z\ = a) = 
and Eqs. ©-(IS]) show that the energy density of the disk is continuously 
matched with the surrounding dust cloud on the vertical borders of the disk. 
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3 Thick disks from the Bonnor solution 



Bonnor jlO] found a simple solution of equations ()2ap -([3]) by starting with 
the following solution of Eq. ([3]) 

2k , , 

^ = Tt^' ^^^^ 

where fc is a constant. The other functions read 



(?k^{r^ + 4^2) 
27rGe'^(r2 + z2)4- 



(19b) 



The energy density is singular at the origin. Now applying the transforma- 
tion z h{z) + b on Eqs. (fT9al) - (fT9bl) and using Eqs. dZ])-®, dSij) and ([TH]) . 
we obtain 

^ = . 2^2 2.rf^\. n2nn { " + + ^r\h + 6)2/j'2}2 
47r2G2e2'^[r2 + (/i + 6)2]io LL v /j v ^ j 



2^2 
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167r2G2e2'^[r2 + (/i + 6) 

'2\r 2 ^(-T , .N2n2 



{(6 + /i)[r2 + (/i + 6)2]/i" 



l-/i'2)[r2-4(/i + 6)2]}% (20a) 



27rGe'^[r2 + (/i + 6)2]4 4vrGe^[r2 + (/i + 6)2 

[r^ - 2(/i + 6)2] [r^ + (/i + 6)^] - (1 - /i'2)[2r^ - 27r\h + hf + 31(/i + 6)^]} 

(20b) 

9c^fc2r2(/t + 6)2(1-/1^2) 

' ~ ' ~ 87rGe-[r2 + (/i + 6)2]5 ' ^^^'^^ 

We nondimensionalise the quantities by defining f = r/a^ z = z/a, h = 
h/a, d = d/a, k = kja^, D = {G^a^/c^)D, S = {Ga^/c'^)S, p = {Ga'^/c^)p, 
k = {Go? /c^)k and Pi = {Go? /c^)Pi. Fig. [T] displays curves of D = as 
function of f and z for the polynomial h with n = 1 and some values of A;, d 
and 6. On the left of each curve the discriminant D is positive. We observe 
that increasing the value of the parameter k keeping the others fixed, the 
level curve D = is moved to the right. The same happens when 6 is 
decreased and the other parameters are kept fixed. 

Fig. [2)[a)-(e) show, respectively, curves of the discriminant Eq. (j20ap . 
the energy density Eq. ()13p . the azimuthal stress Eq. ()14|) . the heat flow 
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function (jlSp and the radial pressure Eq. (|16p . The parameters used are 
n = 1, k = 0.5, d = and 6 = 1. Each curve is a cut taken on the planes 
z = (the curves with largest amplitudes), z = 0.2, 0.4, 0.6, 0.8, 1. From Fig. 
W[i>) we have an annular region of the disk with negative energy density. For 
the azimuthal stress (Fig. [U^c)) there are two rings with pressure (positive 
stress), however near the vertical borders of the disk we only have tension. 
Most of the disk radial extension has heat flow (Fig. (SJd)) that is zero on 
z = 1. From Fig. E^e) we see that the radial pressure (as well as tension 
along z) is concentrated in the central regions of the disk. 



4 Thick disks from a "monopole" and "quadrupole" 
solution 

We consider another solution of Eq. ([3]), namely, a sum of a monopole-like 
and a quadrupole-like term 

V^^^T^ (r2 + z2)5/2 ' ^ ) 
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Figure 2: (a) The discriminant Eq. (|20ap . (b) the energy density Eq. (jl3p . 
(c) the azimuthal stress Eq. (I14p . (d) the heat flow function (llSp and (e) the 
radial pressure Eq. (fT6]) as functions of r for z = 0, 0.2, 0.4, 0.6, 0.8, 1 and 
parameters n = 1, k = 0.5, d = and 6 = 1. In (a) the small plot shows the 
curves for z = 0.8 and z = 1. 
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where ki and k2 are constants (note that ki = 2k when compared with 
Bonnor's solution). Eqs. (|2b[l - (f2c]l and the relation = r^ r give 



N = -p^^ [kiir' + z^f + 2,k,{r^ - Az^)] , (22a) 



r2 



[2fc^(r2 + ^2)4(^2 _ g^2^ ^ 2Akik2{r'^ + z'^fx 



16(r2 +z2)8 

_ 18^2^2 _^ -Lg^4^) _^ 9/^2(^9^6 _ 288r^z2 + 672r2z^ - 256z^)] , (22b) 



,2 



[A:?(r2 + z2)4(r2 + Az^) + 6A:iA;2(r2 + z' 



2\2, 



8^Ge'^(r2 + ^2)8 

(Sr'' + 12r2z2 - lez"^) + 9A;|(9r6 + 72r^z2 - ASr'^z^ + 64^^)] . (22c) 

Applying again the transformation z h(z) + b on the above solutions, we 
obtain for the physical quantities of the disk exact but long expressions that 
are not very illuminating. In Fig. [3] we show how the Eq. D = is changed 
with the parameter k2 = k2/a^ (all other quantities were adimensionalised 
as in Sec. [3]), for /i with n = 1, (i = 0, fci = 1 and 6 = 1. Values for k2 
somewhat greater then 0.053 give rise to a second level curve where Z? = 0. 

Fig. IHa)-(e) display, respectively, curves of the discriminant D the en- 
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Figure 4: (a) The discriminant D (b) the energy density Eq. ()13p . (c) the 
azimuthal stress Eq. ()14p . (d) the heat flow function (jlSp and (e) the radial 
pressure Eq. (fT6]) as functions of f for i = 0, 0.2, 0.4, 0.6, 0.8, 1 and parame- 
ters n = 1, ki = 1, d = and 6=1. In (a) the small plot shows the curves 
for z = 0.8 and z = 1. 
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ergy density Eq. (jl3p . the azimuthal stress Eq. (jl4p . the heat flow function 
(jlSp and the radial pressure Eq. (jl6p . The parameters used are n = 1, 
ki = 1, ^2 = 0.05, d = and 6 = 1. As in Fig. [2]each curve is a cut taken on 
the planes z = (the curves with largest amplitudes), z = 0.2, 0.4, 0.6, 0.8, 1. 
When compared with Fig. [5](a)-(e), the curves are shifted to the right. We 
also note the appearance of a new local maximum point in the energy den- 
sity (Fig. m^b)) and a new local minimum in the azimuthal stress (Fig. IHc)). 
The curves for the heat flow function and radial pressure are not changed 
qualitatively. 



5 Thin disks and the Cooperstock and Tieu galac- 
tic model 

The thin disk limit is achieved by choosing the function h{z) = \z\ in the 
transformation z — > h{z) + b. Since h' = 29{z) — 1 and h" = 2(5(z), where 
9{z) and 6{z) are, respectively, the Heaviside function and the Dirac distri- 
bution, the components (j5ap ~ (f5el) of the energy-momentum tensor will have 
a distributional part with support on the plane z = 0: T^^ = Q'^5{z). The 
nonzero components of are 

Qt = -^,h], Q^ = -T-7t^ 1 + > (23a) 



Qi = 337^1237' = -ZZPTTT. ( + ^,hU (23b) 



S-KGr'^e"' 8^Ge^ V 

The calculation of the eigenvalues result in 



Because the discriminant is always negative, the surface energy density a of 
the thin disk is given \>y a = Sj (2c^) and the azimuthal stress by P^p = —S/2. 
Thus the disk is composed of matter with an equation of state P^p = —c^a. 
If o" > we have tensions and it may be interpreted as an equation of state 
of matter formed by concentric loops of cosmic strings [41]. If o" < we 
have pressure and matter with negative energy density. Objects with this 
equation of state are known in the literature as struts and they appear to 
stabilize certain superpositions of static isolated bodies in General Relativity 
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(see, for instance, In either case the thin disks are made of exotic 

matter with unusual properties. Note that this result is general in the sense 
that it is independent of particular forms for the metric functions and v. 

In the galactic model presented by Cooperstock and Tieu |34| the rele- 
vant field equations that describe the distribution of galactic matter are, to 
order G, given by 

N^rr + iV,.. - ^ = 0, (26) 

r 

P = ^ i^l + Nl) . (27) 

The metric function N{r, z) is calculated from the measured tangential ve- 
locity V = cN/r = c^^r; and then the density distribution is determined by 
the nonlinear relation (I27p . A convenient set of basis functions that are a 
solution of equation (j26p is given by 

N = -Y, Cnkne-'^-^'^rJiiKr), (28) 

n 

where Ji is the Bessel function of order 1 and the C„ are constants. The ab- 
solute value of z must be used to ensure the reflexive symmetry with respect 
to the plane z = 0. However, this is equivalent to make the transformation 
z — > |z| and by the results stated above, introduces on the plane z = an 
additional rotating disk made of exotic matter. It is worth to note that 
in another context de Araiijo and Wang |33] have considered a solution of 
equation ()26p in the form e^'^'r Ji(r) and also comment that this introduces 
an additional mass layer with an equation of state P,p = —(?a and with heat 
flow. 

In principle the exact solution (j28p may also be used to generate thick 
disks via the transformation z h{z) (with 6 = 0), however the calculation 
of the other metric function v using relations (j2bp is far from trivial. But at 
least the zeros of the discriminant equation ([8]) can be determined straight- 
forwardly, since they do not depend on the function u. As an example we 
take as 

10 

AT = -a ^ CnKe-'^'^-'^^'^rJiiaKf), (29) 

n=l 

where the C„A;„ and were fltted for the Milky Way as given in Table 1 of 
|34j . Fig. [5] shows some curves of D = for h with n = 1 and n = 3. The 
variables were adimensionalised as in Sec. [3] and the half-thickness of the 
Milky Way was set to a = 0.75 Kpc. On the left of the curves we have small 
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0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02 

r 

Figure 5: Curves of = as function of r and z for h with n = 1 and 
n = 3. The curves were calculated with Eq. (I29|) fitted for the Milky Way 
(table 1 of tMj). 

intervals of radii where D > 0. In the limit a ^ these intervals would 
shrink to zero, leaving a singular layer on z = with exotic matter. 

6 Discussion 

Using the van Stockum class of metrics that describe spacetimes with sta- 
tionary axially symmetric distributions of dust, and applying the "displace, 
cut, fill and reflect" method, we constructed thick rotating disks surrounded 
by rotating dust. In general the disks have equal radial pressures and ten- 
sions in the z direction. The solution found by Bonnor [40j and a sum of a 
"monopole" and "quadrupole" solution were used to study particular mod- 
els of rotating disks. They were found to have annular regions with negative 
energy density and azimuthal pressure. Heat flow is concentrated near the 
plane z = and is present in most of the radial extension of the disk. These 
particular examples suggest that the "displace, cut, fill and reflect" method 
used on the van Stockum class of metrics will generate disks with some 
unrealistic physical properties and thus are hardly useful in galactic mod- 
elling. However, the aplication of the above mentioned method to solutions 
obtained using a different approach (see, for instance, [38]) may result in 
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more realistic and useful models. 

In the thin limit we obtain disks with heat flow everywhere and com- 
posed of exotic matter with an equation of state P^p = — c^o", which may 
be interpreted as an equation of state of cosmic strings or struts. We also 
comment the new galactic model proposed by Cooperstock and Tieu [M] 
and the presence of an additional singular layer of exotic matter in their 
model. 
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